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We summarise the latest experimental QCD studies based on data from LEP. Measurements 
of the quark and gluon jet fragmentation functions are discussed, including a new algorithm to 
infer the properties of unbiased gluon jets. We describe a new test for destructive interference 
in the radiation of soft gluons from a three-parton system. Finally, we report the latest 
combined value of the strong coupling, measured using event shape observables. 



1 Introduction 

The LEP collider at CERN was used to study e"'"e~ annihilation at centre-of-niass energies in the 
range y/s = 91-209 GeV, during the years 1989-2000. Four multi-purpose detectors (ALEPH, 
DELPHI, L3 and OPAL) were positioned at 90° intervals around the circular accelerator. Events 
of the type e^e~ — > TP /■^ — > hadrons were used extensively to test QCD predictions, and to 
measure the colour factors and strong coupling. In this article, we outline four of the most recent 
studies performed by the LEP collaborations. 

The text is organised as follows: in Sections 2 and 3, we report measurements of quark and 
gluon jet properties, published by the OPAL Collaboration. Section 2 deals with the scaling 
of fragmentation functions over a range of energies, for jet samples with a variety of flavour 
compositions, while Section 3 focuses on a new study of 'unbiased' gluon jets. In Section 4, 
we discuss an innovative test for the presence of destructive interference in the radiation of 
particles from a three-jet system, by the DELPHI Collaboration. Finally, in Section 5, we 
report the progress of a combined measurement of the strong coupling derived from event shape 
observables, performed by the LEP QCD Working Group. 

We emphasise that the results discussed here do not constitute an exhaustive list of the recent 
and ongoing QCD studies by the LEP collaborations. In particular, the results of pentaquark 
searches, colour reconnection studies, and 77 physics will be omitted due to time constraints. 



2 Scaling violations of quark and gluon jet fragmentation functions 

The fragmentation function, D^{x,Q'^), is defined as the probabihty that a parton a, which is 
produced at a short distance of order 1/Q, fragments into a hadron h carrying a fraction x of 
the momentum of a. QCD predicts that the multiphcity of soft gluon emission from a gluon 
source should be higher than that from a quark source, due to the inequality of the colour factors 
Ca and Cp. We therefore expect softer fragmentation functions for gluon jets. Furthermore, 
one can predict the dependences of the fragmentation functions on the scale Q^, by means of 
the splitting functions i'q-^qg ~ Cp and i'g-^gg ~ Ca'- the scaling violations for gluon jets are 
expected to be larger than those for quark jets. 

In a recent study by the OPAL Collaboration,^ the fragmentation functions have been mea- 
sured for a variety of quark and gluon jet samples in e+e~ -^ {TP /^) -^ qq(g) interactions 
at \/s = 91.2 and 183-209 GeV. There are two experimental approaches to the identification 
of jets. Firstly, one can use a jet-finding algorithm (the Durham algorithm in this case): jets 
obtained by this method are biased, because their properties depend on the choice of jet finder, 
and on its associated parameters. Alternatively, jets may be defined as inclusive hemispheres 
of particles in a back-to-back qq or gg system: these unbiased jets correspond to the definitions 
commonly used in theoretical calculations. The new OPAL measurements include seven types 
of fragmentation functions: the udsc, b, gluon and flavour-inclusive fragmentation functions for 
biased jets, and the udsc, b and flavour-inclusive quark fragmentation functions for unbiased 
jets. A previous study of unbiased gluon jets has been performed with rare OPAL events of 
the type e"'"e~ -^ qqgjjjci, in which the gluon 'gind' is identified as the hemisphere opposite to 
a hemisphere containing two tagged b quark jets which are almost collinear.^ The new results 
presented here for biased gluon jets will be compared with the older results for unbiased jets. 
An alternative approach to unbiased gluon jets will be discussed in the next section. 

Three methods have been used to distinguish between quark and gluon jets, and between 
difl'erent flavours of quark jets. In the b-tag method, a neural network based algorithm is used 
to select b quark jets in three-jet events. The untagged jet in an event containing two b-tagged 
jets (or the lowest-energy untagged jet in an event with one b-tag) is selected as a biased gluon 
jet. In events with no b-tagged jets, all three jets are selected as biased udsc quark jets; a 
correction is applied to remove gluon jets from the sample. In the energy- ordering method, the 
jets 1, 2 and 3 of a three-jet event are ordered such at Ei > E2 > E^. Jet number 2 is assigned 
to the biased flavour-inclusive quark jet sample, while jet number 3 is assigned to the biased 
gluon jet sample. Finally, in the hemisphere method, a b-tagging algorithm is applied to an 
inclusive sample of hadronic events. Each event contains two hemispheres, which are regarded 
as unbiased udsc or b quark jets. In all three methods, an unfolding procedure is employed to 
correct for impurities in the selection. 

For experimental purposes, we define the scale Q for a biased jet of energy E^^t to be 
Qjet = -Ejct sin(6'/2), where 6 is the angle to the closest jet. In the case of unbiased jets, we take 
Q = y/s/2. The momentum fraction x is replaced by the quantity Xe = £'h/-E'jct, where -E^ is 
the energy of the hadron and -Ejet is the energy of the jet to which it is assigned. 

Next-to-leading order QCD predictions for the fragmentation functions have been calculated 
by three groups: Kniehl, Kramer and Potter (KKP),^ Kretzer (Kr),'^ and Bourhis, Fontannaz, 
Guillet and Werlen (BFGW).^ The predictions correspond to unbiased jets, and are derived 
using the DGLAP evolution equations, from a set of measured fragmentation functions at a 
fixed input scale. 

The measured scale dependences for udsc, b and gluon fragmentation functions are shown 
in Figures 1 and 2, for different ranges of x^- For the udsc quark fragmentation functions, 
good agreement is found with the NLO predictions, except in the lowest and highest regions 
of Xe- The agreement is poorer in the case of b quark and gluon jets; however, the gluon jet 
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Figure 1: Scale dependence of the udsc and b jet fragmentation functions in different bins of Xe- The 'scale' 

denotes Qjot for the biased jets and v^/2 for the unbiased jets. The inner and outer error bars indicate statistical 

and total uncertainties respectively. The data are compared to NLO predictions by KKPf Kr '' and BFGW.^ 
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Figure 2: Scale dependence of the gluon jet fragmentation functions in different bins of a::E- The 'scale' denotes 

Qjet for the biased jets and Ejct for the previously published results using unbiased 'ginci' jets. The results for 

biased jets are shown separately for the the b-tag (BT) and energy-ordering (EO) selection methods. The inner 

and outer error bars indicate statistical and total uncertainties respectively. 

fragmentation functions do exhibit stronger scaling violations than the quark jets, as expected. 
Good agreement is found between the biased and unbiased jet samples, suggesting that Qjet is 
an appropriate scale in events with a three-jet topology. The results obtained using the b-tag 
and energy-ordering methods are also consistent with one another, and with previous results 
from DELPHI and OPAL. 



3 Studies of unbiased gluon jets using the jet boost algorithm 



As we have mentioned in Section 2, a distinction exists between the biased jets obtained from 
experimental data using a jet-finding algorithm, and the unbiased jets used in theoretical calcu- 
lations. An unbiased gluon jet would correspond to one hemisphere of a back-to-back gg system, 
which is not seen in e'''e~ annihilation. Instead, rare events of the type e^e~ -^ qqginci have 
been used, in which the q and q jets are approximately collinear, leaving an unbiased 'gmci' jet 
in the opposite hemisphere.'^ Another OPAL study has used a more indirect method, whereby 
the results obtained from two-jet qq events are subtracted from those in qqg events.^ Unbiased 
gluon jets have also been obtained from radiative T -^ 7gg decays at CLEO.'^ Recently, however, 
a new approach known as the jet boost algorithm has been considered in e^e~ annihilation.^ In 
this method, a qqg system is decomposed into two independent qg and qg dipoles. The dipoles 
are boosted into a symmetric frame, such that the angle 2a between the q and g is the same 
as that between the q and g, as shown in Figure 3(a). Further Lorentz boosts (3 = cos a are 
then applied independently to the two dipoles, such that they are each back-to-back, as in Fig- 
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Figure 3: The jet boost algorithm 
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Figure 4: The mean charged particle multiplicity of unbiased gluon jets, (rig^'uon) i ^^ ^ function of the jet energy E* 



ure 3(b). Finally the two dipoles in their different frames can be recombined, yielding an event 
with the colour-structure of a gg system in a colour singlet state, as shown in Figure 3(c). The 
hemisphere containing the gluon in this event corresponds to the theoretical definition of an 
unbiased gluon jet. Unlike the unbiased quark jets discussed in Section 2, the energies E* of 
these gluon jets are not fixed by the e'''e~ centre-of-mass energy. 

The jet boost algorithm has recently been studied by the OPAL Collaboration, using e'^e" -^ 
qqg events at the Z^ resonance.^ Gluon jets are selected using a combination of b-tagging and 
energy-ordering; cuts then are imposed on the softness and collinearity of the quark and gluon 
jets. Using the HERWIG Monte Carlo event generator, it has been established that the unbiased 
gluon jets obtained using the jet boost algorithm should have properties consistent with a true 
back-to-back gg pair. It has also been shown that the properties are essentially independent of 
the jet-finder used to identify the q, q and g jets, as expected for an unbiased jet. 

The mean charged multiplicity of the unbiased gluon jets, {n'^^^^), is shown in Figure 4, as 
a function of the jet energy E* The results are found to be consistent with previous OPAL 
measurements based on the 'gincl' and 'subtraction' methods,^'^ and with the HERWIG prediction 
for genuine gg pairs. QCD evolution fits are also performed for two calculations: one based 
on next-to-next-to-next-to- leading order (3NL0) perturbation theory for a scale-dependent as, 
and the other based on an exact solution for a fixed as- Two free parameters are used in 
each case, and no hadronisation corrections are applied. The QCD scale parameter A is hence 
found to be A = 0.296 ± 0.037 (total) GeV, compared to the corresponding quark jet result 
A = 0.190 lb 0.032 (stat.) GeV.^ Similar analyses have been performed for the first two non- 
trivial factorial moments of the charged particle multiplicity distribution, -F2, gluon and i*3,giuon? 

Fragmentation functions have also been obtained for the unbiased gluon jets, at energies 



£'*=14.24 and 17.72 GeV. Since the OPAL Collaboration has previously measured the fragmen- 
tation functions for unbiased gluon jets at 40.1 GeVf and for unbiased quark jets at 45.6 GeV,^° 
the DGLAP evolution equations can be used to construct a QCD prediction for the new measure- 
ments at lower energies. Such predictions are dependent on the strong coupling Og, which can 
therefore be extracted from a one-parameter fit. The fits are in good agreement with the mea- 
sured fragmentation functions, yielding as(Mz) = 0.128 it 0.008 (stat.) it 0.015 (syst.). Although 
not competitive with other measurements of the strong coupling, this result is compatible with 
the world average, and has provided a unique consistency test of QCD. 

4 Coherent soft particle production in Z° decays into three jets 

Interference effects are fundamental to all quantum mechanical theories, including the gauge 
theories of the Standard Model. Evidence for coherence effects in QCD comes, for example, 
from the "hump-backed plateau" in the logarithmic scaled momentum spectrum of hadrons, due 
to suppression of low energy particle production, and from the string effect in three-jet e^e~ 
annihilation events.^'^ However, there are arguments against the conclusiveness of this evidence.^^ 
In a new study by the DELPHI Collaboration,^^ a direct test is made for the presence of 
destructive interference in the radiation of soft gluons from a three-jet system. When a quark 
emits a hard gluon at a small opening angle, the quark-gluon system may behave on large 
distance scales as a single entity. A soft gluon radiated perpendicular to the q~g pair will not 
be able to resolve the individual colour charges of the quark and gluon; it should therefore 
be regarded a coherent emission from the parton ensemble. The leading order cross section, 
d(T3, has been calculated for coherent soft gluon emissions perpendicular to the plane of a qqg 
system: ^^ 



d(T3 = d(T2 



(1) 
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ncr -I- ncr — - 

ACf 

Here d(T2 represents the corresponding cross section for soft gluon emissions perpendicular to 
the axis of a two-jet (qq) event, and the antenna function ij for a pair of partons (i,j) with 
an opening angle Oij is defined by ij = 1 — cos 6ij. The last term of Equation (1), inversely 
proportional to the square of the number of colours A'c = 3, is due to destructive interference. 
By measuring the ratio d(T3/d(T2 in samples of two- and three-jet events, one should be able to 
verify this (l/N^) qq interference term, and to determine the colour factor ratio Ca/Cf based 
on the leading order expression given above. 

Hadronic Z^ decays were recorded by the DELPHI detector in e"''e~ annihilation at ^/s = 
91 GeV. The angular ordered Durham algorithm ^^ was used, with a fixed resolution parameter 
ycut = 0.015, to determine the number of jets in each event. For each three-jet event, the charged 
particle multiplicity N^ was measured in a 30° cone oriented perpendicular to both sides of the 
event plane. The corresponding multiplicity A'^2 was also measured in two-jet events, for a cone 
perpendicular to the event axis; the azimuthal angle of the cone was chosen randomly in this 
case. In Figure 5, we show the three-jet cone multiplicity N^ as a function of the inter-jet 
angles 62 and ^3, where 9i is the angle between the two jets opposite to jet i, and the jets are 
energy-ordered such that Ei < E2 < E^. From the measured value of N2, a prediction can be 
calculated for iV3 using Equation (1): these predictions are shown in Figure 5, with and without 
the destructive interference term. The expression without the (1/A'^)qq term is incompatible 
with the data, while the fully coherent prediction is in good agreement. If the ratio of colour 
factors, Ca/Cf, is fitted to the data, one obtains the value 2.211 it 0.014 (stat.) it 0.053 (syst.) 
with x^/ndf = 1.3. Although this result is valid only at leading order, and does not include an 
estimate of the theoretical uncertainty, it is in astonishingly good agreement with the expectation 
Ca/Cf = 2.25. The leading order QCD prediction for the soft gluon multiplicity, including 
destructive interference, has therefore been verified convincingly by the data. 
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Figure 5: Particle multiplicity in cones of 30° opening angle perpendicular to the three-jet event plane, as functions 
of the opening angles 62, 0^. The inner error bars are statistical, while the outer bars include both statistical and 

systematic uncertainties. 



5 A combined measurement of as using event shape observables 

Event shape observables have been used extensively to test QCD predictions, and to measure the 
strong coupling as, in both e''~e~ annihilation and deep inelastic scattering.^^ The LEP collabo- 
rations have published measurements of as for events in the energy range y/s = 91-209 GeV,^^ 
using the following six event shape observables: thrust (T), heavy jet mass (Mh), C-parameter, 
total jet broadening {B^), wide jet broadening (i?w)) ^-iid the Durham two-to-three jet res- 
olution parameter (^3). For each distribution, an 0{a^) perturbative prediction is matched 
with an NLLA resummation and corrected for hadronisation effects using the PYTHIA Monte 
Carlo event generator. The QCD predictions are then fitted to the LEP data, with as as a 
free parameter. The LEP QCD Working Group has combined these into a single preliminary" 
measurement of as at the TP mass scale, with particular attention to estimation of theoretical 
uncertainties ^^ and to the treatment of correlations.^'^ The result is 

as(Mz) = 0.1202 ± 0.0003 (stat.) ± 0.0009 (expt.) ± 0.0013 (hadr.) ± 0.0047 (theo.) , 

where the three systematic uncertainties are due to experimental effects ('expt.'), Monte Carlo 
hadronisation corrections ('hadr.') and higher-order terms in the perturbative QCD predic- 
tions ('theo.'). Our measurement is in good agreement with the current world average.^^ 

°A publication is expected later in 2004, when measurements provided by the collaborations have become final. 



"s((3) 

0.125 
0.120 
0.115 
0.110 
0.105 
0.100 



<- as{Q) fit at single energy 



<^ as(Q) global fit 

'^ ± 0"stat. 




LEP Preliminary 



100 



150 



200 

Q (GeV) 



1-T 



Mg/,s- 



C 



St 



-Bw I ^*^- 



- My 3, 



All 



_L 



_L 



-l-»H 1 



_L 



_L 



0.110 0.115 0.120 0.125 0.130 

as {Mz 



Figure 6: Left figure: Measurements of the strong coupling as{Q) using event shape observables at LEP. Each point 
is a combination of all available measurements at a single energy; the inner error bar represents the combined 
statistical uncertainty, while the outer bar is the total uncertainty. The curve indicates the running of as{Q) 
predicted by QCD at three- loop order, based on our combined measurement of as at the Z" mass scale, as(Mz). 
The statistical and total uncertainties in as{Mz) are indicated by the dotted curve and grey band respectively. 
Right figure: Combined measurements of Qs(Mz), using individual event shape observables. 
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